Commensurate and incommensurate correlations in Haldane gap antiferromagnets 
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We analyze the onset of incommensurabilities around the VBS point of the S — 1 bilinear- 
biquadratic model. We propose a simple effective field theory which is capable of reproducing 
all known properties of the commensurate-incommensurate transition at the disorder point v bs- 
Moreover, the theory predicts another special point #dis P , distinct from the VBS point, where the 
Haldane gap behaves singularly. The ground state energy density is an analytic function of the 
model parameters everywhere, thus we do not have phase transitions in the conventional sense. 
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I. INTRODUCTION 



The S = 1 bilinear-biquadratic chain 

H = [cosflS^+x + sm6 (SA+i) 2 ] 



(1) 



is one of the prototype models for the physics of Haldane 
gapd antiferromagnets. Its zero temperature phase dia- 
gram phas been the subject of intensive studies in recent 
yearsJa By now it is well-established that the energy 
gap (Haldane gap) persists in a wide range —it/ 4 < 8 < 
7r/4 around the conventional Heisenberg point 6 = 0. 
The model is gapless at the special pointsa 8 = ±7r/4 
beyond which other phases with qualitatively different 
physical properties appear. Although the gap and the 
hidden (string) ordera characterizing the Haldane-gap 
systems persist for the whole Haldane phase — ir/4 < 8 < 
7r/4, one can divide this interval into (at least) two, some- 
what different subphases. ul] These subphases are sepp. 
rated by the so called valence-bond-solid (VBS) pointE£l 
#vbs = tan -1 1/3 w 0.10247T where the ground state prop- 
erties can be obtained exactly. The two subphases differ 
in the form of the long distance asymptotics of the two- 
point correlation function (SfSf: n ). In the "commensu- 
rate" Haldane phase (C-phase) for — 7r/4 < 6 < # v bs the 
leading behavior is expected to be 



(s?s? +n ) ~ (-iy 



,1/2 



for n 



oo, 



(2) 



while in the ' 
phase) for 6> v b s 
the form 



incommensurate" Haldane phase (IC- 
< 8 < 7r/4 this was predicted to take 



cos(qn 



for 



(3) 



where q = q{8) G (tt, 2ir/3) is a ^-dependent incommen- 
surate wavenumber, is a phase shift, and £ = £(0) is the 
correlation length. At the "C-IC transition point" (also 
called the "disorder point") 8 clc = # v bs the correlation 
function is known rigorously and it is purely exponential 
without any algebraic prefactor 



l/ < - fa, 



(4) 



with ^(^vbs) = l/ln.3 ~ 0.9102. Correlation functions of 
some other operators can also be studied rigorously at the 
VBS pointE In particular, ((S z ) 2 (S z ) 2 +n ) has a similar 
purely exponential decay with £=l/ln3at# = # v bs- 

The commensurate-incommensurate (C-IC) transition 
of the spin-1 chain bears much similarity to C-IC transi- 
tions found in other models, e.g., in the anisotropic 2D 
Ising model on the triangular lattice at finite tempera- 
ture, wheae it can be analyzed rigorously using the exact 
solutionis Other, higher dimensional, not exactly solv- 
able models with disorder points were investigated using 
an RPA approach to the susceptibility in Ref. [l^. In gen- 
eral, C-IC transitions can be divided into two categories 
(two kinds). In this classification scheme the transition 
at the VBS point of the S — 1 chain is a C-IC transi- 
tion of the first kind, with the property that the incom- 
mensurate wavenumber q in the IC regime is parameter 
dependent. For C-IC transitions of the first kind the cor- 
relation length is predicted to behave on the C and IC 
sides of the disorder point 8 c i c as 



(16 



= finite 



(5) 



ic 



with £(0 C ic) 7^ at the transition point. The characteris- 
tic wavenumber is expected to vary on the IC side as 

- g(0cic) ~ \d-0 cic \ 1/2 . (6) 

Moreover, the RPA theoryQ also predicts the change of 
the asymptotic form of the correlation functions at the 
disorder point: the algebraic prefactor is n- {D - 1)/2 in D 
dimensions except at the disorder point where D — > D' , 
D' < D should be taken, reflecting a "dimensional re- 
duction". In our case D = 1 + 1 = 2, £)' = 1 as is 
seen in Eqs. (§-|]). These general features of the C-IC 
transitions of the first kind have been tested numerically 
for the spin-1 chain, and except for some observed devi- 
ation from the predicted form in Eq. (^) slightly above 
6* v bs, all were justifiedJa Note that around the VBS point 
finite size corrections are very small, thus the numerical 



1 



results (exact diagonalization and DMRG) are extremely 
precise. 

Within the IC subphase some other special points 
can be defined. The first one is 6>di sp ~ 0.12107T, 
where the second derivative of the magnon dispersion 
at k = 7r vanishes, and the dispersion becomes quartic. 
When the chain is subject to a uniform magnetic field 
the magnetization-vs-field curve has an aapmalous, non- 
square-root-like singularity at this point £3 For 9 > 9<xi S p, 
the second derivative at k = tt is negative and the 
gap takes its miftimum value for a momentum differ- 
ent from k — ttB In this region a strong enough field 
causes the Haldane gap to collapse into a iwo-component 
Luttingcr liquid (LL) phase instead of a conventional 
one-component LL phase.E3 The next special point is 
$max ~ 0.1237T where the Haldane gap takes its maximum 
value. Note that this point has less physical relevance, 
since it strongly depends on the actual definition of the 
Hamiltonian; a ^-dependent rescaling of the model can 
easily change its location. Finally one can define the Lif- 
sitz point $Lif s ~ 0.13147T, where the incommensurability 
manifests itself in the structure factor. This is different 
from the disordered point iu.massive models due to the fi- 
nite linewidth of the peaks.U Figure [j] shows these special 
points. 



play a role at low energies. Many propertieSp.oi-.the 
system, from ground-state correlation functioiplpla to 
the onset of magnetization in uniform fieldsr 3 H 19 H 2C i can 
be extremely well approximated using a massive rela- 
tivistic free boson theory. Such an approximate the- 
ory can be derived directly from the non-linear er-model 
(NLctMj) description of the spin chain at the Heisenberg 
pointjlfl or from the Majorana fermion representation of 
the integrable_Takhtajan-Babujian model in the vicinity 
of 6 = — 7t/4.Ej However, all these microscopic theories 
utilize the a priori assumption that the important low- 
energy fluctuations are at momenta k — and tt, and 
thus they cannot account for the C-IC transition, nor 
can they give any reliable description of the IC regime. 

Although we do not have a rigorous microscopic justi- 
fication all results we possess are consistent with the as- 
sumption that the elementary excitations are essentially 
free bosons in the IC regime, too. Thus the aim of the 
present paper is to extend the free boson description to 
the whole Haldane phase, and give a general theory which 
is capable of accounting for the C-IC transition in simple 
terms. In the lack of a detailed microscopic formulation, 
however, this theory, at present, is only phenomenologi- 
cal. 



1 -component 
Luttinger liquid 



2-component 
Luttinger liquid 




C - Haldane 



IC - Haldane 



e 

FIG. 1. Schematic phase diagram showing the neighbor- 
hood of the C-IC transition with the special points. The crit- 
ical magnetic field, where the Haldane gap collapses, equals 
the value of the energy gap. 

In both the C and IC Haldane phases the energy spec- 
trum consists of a discrete triplet branch of "magnons" 
separated by finite gaps from the singlet ground state, 
and also from a continuum of higher lying "multi- 
magnon" excitations in a wide momentum range. This 
one-magnon branch is clearly discernible around the 
edges of the Brillouin zone k = ±7r. However, it merges 
into the continuum and vanishes due to magnon scat- 
tering processes in an extended range around k = 0. 
The energy spectrum, especially in the Heisenberg point 
and the n |iBS. jQoint was studied numerically by many 
authors .Q'EjO'Ej All concluded that the magnon-magnon 
interactions are rather weak, and bound states do not 



II. EFFECTIVE THEORY - CONTINUUM 
VERSION 

A continuum field theory to describe Heisenberg anti- 
ferromagnetic chains (8 = 0) with integer S was devel- 
oped by Haldaneill This is the nonlinear cr-model (NLcrM) 
- without topological terms - which can be derived in the 
large S limit, but whose implications are believed to hold, 
for S = 1 too. The NLctM is defined by the LagrangianEJ 



C 



1 



-(5 t 0) 2 - v(d x 4>f 



(7) 



where 4>(x,t) is a vector field with unit length (jr = 1, v 
is a nonuniversal constant (velocity) setting the energy 
scale, and g = 2/S is the coupling constant. The associ- 
ated Hamiltonian is 



v 

H =2 



9 



1 



(8) 



where the momentum canonically conjugate to (p turns 
out to be 



l(x,t) = — 
gv 



(f> x 



dt 



(x,t). 



(9) 



The spin operator S n can be expressed in terms of <f> and 
I as 



S n {t) = {-l) n S$(x n ,t) + 8xl(x n ,t) 



(10) 



where x n = nSx, Sx being the lattice constant (usually 
set to unity), and S = 1 in the present case. At the 
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mean field level the NLtrM can be well approximated by 
a massive, essentially free vector-boson theory 



1 
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c = -{d t 4>y - v{d x (j)Y - m d (t> 

v 



(11) 



now without any constraint on the </> field. The boson 
field <j> varies smoothly on the scale of the lattice con- 
stant in the commensurate regime (the pure Heisenberg 
model, which the theory applies for, is in the C phase) 
and thus higher derivatives, neglected in Eq. (|ll|), are 
indeed small. The Lagrangian gives rise to a relativistic 
dispersion cu(k) — y/ m 2 + v 2 k 2 /2m which takes its min- 
imum at k — 0. [Note the factor (—1)™ introduced in Eq. 
( |l0| ) shifting all momenta by it.] 

The NLcM summarized above gives a good descrip- 
tion of the low-energy behavior of the S = 1 Heisenberg 
chain, but is unable to describe the C-IC transition in 
the bilinear-biquadratic chain. The key feature missing 
is that in the vicinity of our special points the shape of 
the magnon dispersion changes drastically, and the min- 
imum at k — 7r (k = in the momentum-shifted boson 
language) splits. Obviously, such an effect can never be 
obtained from a relativistic field theory such as Eq. (|ll| ) 
and we need to consider a non-relativistic model. The 
simplest continuum Lagrangian with this property is 



C = (d v 



2 a2 



(12) 



At this level a, b, and m are (9-dependent phcnomcnolog- 
ical parameters. For stability b is supposed to be posi- 
tive, and a is assumed to change sign at the special point 
#dis P - [Note that in the relativistic case b = 0, and a is in 
fact v 2 , where v is the "speed of light". Formally taking 
a < 0, which we will consider in the sequel, implies an 
imaginary v. This only means, however, that one should 
redefine the physical (quasi)particles to live around the 
two new minima of the dispersion. In fact, one should in- 
troduce two new particles, now with v =real, for the two 
minima. This immediately leads to a two-particle (two- 
band) description which was analyzed in detail in Ref. 
for the massless case (above the critical magnetic field 
and was shown to result a two-component (two-band) 
Luttinger liquid there. Here we shortcut these complica- 
tions by simply assuming that the two particles are two 
chiral components (k < and k > 0) of the same boson 
with a < 0.] 

Since our phenomenological model is not derived di- 
rectly from the microscopic Hamiltonian in Eq. (Q) the 
connection between the boson field and the original spin 
variables is not known rigorously. In principle S n can be 
expanded in powers of the field (j) and its derivatives, and 
any term permitted by the symmetries of the model can 
appear. Being guided by the NLcrM description at the 
Heisenberg point, in the following we will assume that 
the first two terms (linear and quadratic in (f>) are 

Suit) = g 4> (-l) n $(x n ,t)+g l l(x n ,t)+0{\$\ 3 ), (13) 

where and g% are unknown constants. Note that the 
I term, defined by Eq. (^), is the most relevant 2-boson 



term which is even under parity and odd under time re- 
versal as it should be. Although, we cannot exclude the 
possibility that up to 0(4> 2 ) some other terms with higher 
derivatives also appear, such extra derivatives, due to 
the finite mass gap, would not modify the long distance 
asymptotics of (SoS n ) . Equation ( p"3| ) shows that in the C 
phase where spin-spin correlations are antiferromagnetic 
(j) varies smoothly, and higher order derivatives in the 
Langrangian have minor role. Around the C-CI transi- 
tion point, and in the IC-phase, however, this is no longer 
the case: as S n picks up incommensuration, 4>(x) must 
do so, as well. It no longer varies smoothly, and higher 
derivatives in the Lagrangian cannot be neglected. This 
is the effect which wc try to take into consideration by 
the b term in Eq. jl^). 

The asymptotic behavior of the correlation function 
(SfS* +n ) is determined by the one-boson term; the two- 
boson term and in general any multi-boson terms only 
constitute minute corrections which decay at least twice 
as rapidly. Thus, up to the two-boson term in Eq. jl^ ) 
the spin-spin correlation function has the behavior 

(S^(t)S^0)) =gl(-) n G4x n ,t) + g?G l (x n ,t), (14) 



with 



G,(x,i) ee (l*(x,f)i*(0,0)>. 



(15) 
(16) 



These are the quantities we will now calculate. 

The Euler-Lagrange equation associated with our La- 
grangian gives a generalized Klein- Gordon equation for 
each component a = x, y, z of 4> 



-ad 2 <j) a +bd: 
whose Green's function is 



4ia 



m 2 4> a 



o, 



G<j,(oj, k) 



— ok 2 — bk 4 — 



IE 



This defines the non-relativistic dispersion 
u(k) = \/m 2 + ak 2 + bk 4 , 



(17) 



(18) 



(19) 



which reduces to the relativistic dispersion of Ref. |l8| 
when 6 = 0. 

The generalized theory can be quantized identically 
to the relativistic Klein-Gordon theoryH The field </> Q , 
a = x,y, z, has the following mode expansion 



4> a (x,t) 



dk 



^i-Kcoik) 



(20) 



where K • X = tot — kx, and the normalization is 
[rff,d?/] = Sp a S(k — k'). Using the mode expansion the 
equal time expressions G^,(x\- and Gi(x) can be easily 
reduced to Fourier transformala 
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G^x) 



Gi(x) 



dk e 



ik:l 



47r u(k) ' 
dk 



(21) 



4tt 



co(k')e 



/!■'.! 



dk e ikx 1 2 
4^w(fc) 2 ^ 



Note that the integral determining G$ also appears as a 
multiplicative factor in G;. 

Since the asymptotic behavior is determined by Ga, 
we start our analysis with this. The evaluation of the 
Fourier transform starts with locating the zeros of ui(k) 
in the complex plane. The four zeros are given by 



k = ± 



— (—a ± \/ a 2 — 4m 2 b 
2b \ 



(22) 



and depend on the single parameter 9, via a, b and m. 
Since our phenomenological model is not derived directly 
from the microscopic one, we have to make several as- 
sumptions about the ^-dependence of a, b and m. We 
will assume that this dependence is smooth (analytical), 
m(9) and b(9) are nonnegative whereas a decreases with 
increasing 9 and changes sign at 9^ sp . We introduce the 
discriminant 



D{9) = a 2 - 4to 2 6, 



(23) 



which is hence an analytic function of 9, and we sup- 
pose that it is positive for 9 < 6* v b s , vanishes at the VBS 
point, and it is negative from 9 v h s to 9 ~ 7r/4 where 
it vanishes again. As we show below, under the above 
hypotheses the phenomenological model provides the ex- 
pected asymptotic behavior of the correlation function in 
all parts of the Haldane phase. 
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— # v bs: 

At this point a > and D = 0. The expression under 
the square-root in ui(k) is a complete square, thus uj be- 
comes quadratic in k and has two purely imaginary zeros, 
±i^2m 2 /a = ±i^Ja/2b [see Fig. |(b)]. The contour of 
integration can be closed in the upper half plane and the 
result is 



G^x) 



1 



2J2a 



(24) 



Because of the purely exponential decay, this indeed cor- 
responds to 9 — 6* v bs- 

< Ovbs'- 

In this region a > and D > 0. We may suppose b > 0; 
the case when 6 = can be obtained by continuity. We 
get four purely imaginary zeros ±iv± where 



v± 



Yb( a± ^. 



Now 



oj(k) = Vl(k 2 + v 2 + )? {k 2 + v 2 _y- 



(25) 



(26) 



is single- valued in the complex plane with two cuts, one 
between — iv+ and — iv- and another one between iv- 
and iv + [see Fig. f§(a)]. (We use the convention that 
z 1 ! 2 has a cut along the negative real axis.) The contour 
of integration can be closed in the upper half-plane and 
drawn onto the upper cut, giving 



G^x) 



2irVb 



-t\x\ 



[{t 2 - v 2 _){v\ -t 2 )\ 



-dt 



(27) 



This function is positive for all x, so with the factor (—1)" 
introduced in Eq. (jl^) we obtain the expected antiferro- 
magnetic modulation. As goes to infinity, the main 
contribution to the integral is coming fcom the vicinity 
of the end point V- , and it is legitimated to expand the 
integrand around this point. For \x\ > |D| -1 / 2 this yields 



Ga(x) 



-v_\x\ 



2 [27 



D * 



Q(\x\ 



(28) 



Together with Eq. (13), we find the expected asymptotic 
form Eq. (||) of the spin-spin correlation function with 
a correlation length £ = l/v_ which is continuous at 

9 = 6*vbs. 
> Qvbs'- 

Here D becomes negative and ui has four complex zeros 
kj = ±u ± iv (J = 1, . . . , 4) with 



FIG. 2. Roots of ui(k) in the complex k plain, (a) 8 < (9 v bs, 
dvhs, (c) 6 > # v bs- Branch cuts are chosen according 



0)6 

to the conventions described in the text. 
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(29) 



If we write oj in the form 



4 



£■>(*) = V^JJ(fc-fc,-)' 



(30) 



we see that it is single- valued on the complex plane with a 
cut between — u — iv and u — iv and another cut between 
— u + iv and it + zt> [see Fig. 0(c)] . For k real we get 
back the original positive function. The integration can 
be carried out along a contour which starts at — u + ioo, 
goes vertically down to — u + iv, passes below the upper 
cut and goes vertically to u + ioo. Next, we replace the 
integral below the cut by an integral going above the cut 
and in the opposite sense, and this latter by the sum of 
the two integrals along the vertical half-lines. These are 
complex conjugate to each other, so finally we obtain 



G^x) 
1 



(31) 



e-^ l cos[u\x\ - if(t)) dt 



irVbJv Vi 2 - v 2 [Au 2 + (t - v) 2 } i [Au 2 + (t + v) 2 } i 
where 



If t — v t + v 

(pit) — — arctan h arctan 

^ w 2 V 2u 2u 



(32) 



Now G$(x) changes sign periodically, and we can identify 
7r — u with the wave number q of the incommensurate os- 
cillation. At the VBS point u = 0, and as it is shown by 
Eq. (|2^), the assumed analyticity of m, a and b assures 
that it has a square-root-type singularity above the VBS 
point in accordance with Eq. (||). The large- 1 x | asvmp- 
totics of G^(x) can be obtained by Watson's lemmaE^I or 
by a direct expansion, 



G^x) = 

e -v\x\ cos ( u |a;| _ I arctan -) 



(33) 



(2n)i 



-mD 



0{\x\ 



Formulas (31) and (|33j) apply for any 9 between 9 V \, S 
and 7r/4, including #disp which is a symmetry point of 
the domain of incommensurate oscillations {u = v). As 
9 approaches 7r/4, D goes to zero and the zeros of u>(k) 
tend to the real axis. Thus, v goes to zero and the cor- 
relation length diverges. This is what we expect at the 
boundary of the Haldane phase 9 — ir/A where the gap 
disappears. At thisjjspecial point the ground state has a 
tripled periodicity,lffl implying u(9 = 7r/4) = 2tt/3. 

For any fixed x, G^,(x) depends analytically on 9 in- 
side the whole Haldane phase. This can be seen from the 
original form Eq. ( pl| ) of G^ (x) by inserting the original, 
non-factorized expression Eq. ([19]) for u>(k). A proof can 
be found in p4| . The argument makes use of the continu- 
ity of the integrand in k real, its (supposed) analyticity 
for any fixed k as a function of 9 in suitable complex do- 
mains, and the uniform convergence of the integral for 9 
in any of these domains. It is interesting to examine an- 
other kind of asymptotics, valid in a close neighborhood 
of v bs, when \D\/a 2 <C 1. In this case Eq. (27) reduces 
to the form 



G$(x) w — / cosh 





x\ sin a 



da 



(34) 



where Io is the zeroth order Bessel function. We can 
arrive at the same equation from Eq. (|3l|), by changing 
the contour of integration (integrating around the upper 
cut). Now analyticity at the VBS point is manifest, be- 
cause in the expansion of the hyperbolic cosine about zero 
only the even powers of \[T) appear. Equation (|34|) shows 
that the crossover to the decay with the \x\~^ prefactor 
sets in at the characteristic distance |a;| ~ D~? : which 
diverges at. the VBS point. This explains the numerical 
difficultieaJ verifying the expected asymptotic behavior 
very close to the VBS point. 

At the VBS point there is an infinite jump in the 
derivative of the correlation length, as predicted by Eq. 
(|). Indeed, Eq. (§f) yields 



d9 



D' 



D-> ~-(9 vhs -9)-i (35) 



3 _o 4y2am 

because D'(# v bs) < 0. On the other hand, from Eq. (|29|) 



dO 



1 / ml a' m? 
h0 2y/a \ 2m 2a a 2 



(36) 



which is finite. The singularity of the correlation length 
at # v bs is in no contradiction with the analyticity of G(x) 
at a fixed x. Indeed, the divergence of the derivative of 
£ was extracted from the single-exponential asymptotic 
form Eq. ( p8|) which, again, is valid only for |i| > (v+ — 

To see the role of the two-boson term G/ in the cor- 
relation functions, we can use the identity (after proper 
regularization) 



dk 

—uj(k)e-' kx 

4-7T 



LU 2 \ I 



d 



dk e- lkx 



dx J J Att u)(k) 



(37) 



where u 2 (id/dx) is a shorthand for m 2 — ad 2 /dx 2 
bd 4 /dx 4 in the present case. With this 



G,l.r) =G,,i.r)^[i^-)G,A.r). 



(38) 



where we have neglected the singular, delta-function term 
of Eq. (|2l]). This term can be evaluated directly for large 
x, knowing the asymptotic form of G^ in the different 
regimes. Using Eqs. (M), ( p8|) and (|33| ) we obtain 



Gi(x) 



if 9 < 9 v hs 
if 9 = 9 v bs 



(39) 



[ci + c 2 cos(2iia; + a)] if 6 > 



vbs i 
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where the constants c\ , ci and a can be expressed 
straightforwardly with a, b and m. It is interesting to 
remark that Gi is exactly zero for any x > at the dis- 
order point, i.e., the 2-boson processes do not contribute 
to the equal time correlation function there. This can 
be easily verified by calculating lo 1 '(id / 'dx)G \p(x) using 
Eq. @ and the fact that D = at the VBS point. Here 
(S*Sq) only contains the term in full accordance with 
the exact solution. For other values of 9, Gi decays twice 
as rapidly as G^. 

It is also of interest to see what predictions our simple 
field theory gives for the analytic properties of the ground 
state energy density and the gap. Using the mode expan- 
sion in Eq. (p0|), the Hamiltonian can be written as 



H 



dk u>(k) 



44 



1 



(40) 



where A is an appropriate UV momentum cutoff, propor- 
tional to the inverse of the lattice constant. From this the 
ground state energy is 



1 f 

E=~ dkuj(k). 

2 J-A 



(41) 



The ground state energy depends analytically on a, b and 
m whenever the zeros of u>(k) are not on the real axis. 
Together with the supposed analyticity of a(8), etc., this 
meansj-khat E(9) is also analytic inside the whole Haldane 
phaserj A straightforward expansion around # v bs yields 



lim [E(6) - E(6 



vbs ) 



(42) 



E 



1 (4ti-4)! 



87ra 3 / 2 



V2b ^(2n-2)!n!(n-l)! 



which is convergent if \D\/a 2 < 1. We notice that in 
general E(ff) can be expressed in a closed form in terms 
of elliptic integrals of the first and second kind. 
The energy gap of the model is by definition 



A = mhifc u>(k). 



(43) 



This is obviously analytic at clc = # v b s but has a singu- 
larity at #disp> where the minima of uj(k) move away from 
k = as the parameter a = changes sign. While for 
9 < #disp the minimum is taken at k = 0, for 9 > # d isp it 
is taken at k = ±y / —a/2b. The gap A and its derivatives 
with respect to 9 on the two sides of # d i sp turn out ot be, 
resp., 



A(# disp - 0) = m, A(9 disp + 0) 



A'(0 disp - 0) = m', A'(# disp + 0) = m', 
A" (9 disp - 0) = m" , A" (9 disp + 0) = m" 



(44) 



46m 



We see that there is a discontinuity in the second deriva- 
tive. This behavior of the effective theory seems consis- 
tent with the numerical results shown in Fig. 9 of Ref. 



III. EFFECTIVE THEORY - LATTICE VERSION 

The effective theory presented above is capable of pro- 
viding a complete qualitative description of the C-IC 
transition of the spin-1 bilinear-biquadratic model in ac- 
cordance with the available numerical data. However, 
in order to give quantitative predictions, too, the the- 
ory needs some refinement. We have seen earlier that 
our simple theory identifies the disorder point 9 cic with 
the point where the discriminant D defined by Eq. ( |2^ ) 
vanishes. It is easy to verify that the second and fourth 
derivatives of u(k) at k = are, resp., 



d 2 uj ,, . a 

aF<* = °> = m' 



d 4 Lo,, , 3D 

w^ k=0) = -^ (45) 



thus in the above theory the fourth derivative vanishes at 
the disorder point. In contrast with this, Golinelli et al.til 
measured numerically the second and fourth derivatives 
at the known disorder point 9 v b s and found 



^(k = 0)- 0.9778(1), Sj(* = 0) 



dk 2 



8k 4 



-1.202(1). 

(46) 



Note that the fourth derivative is only zero far inside 
the IC regime, which seems inconsistent with the above 
theory. 

One step to improve the theory is to realize that the 
model is defined on a lattice, and thus the dispersion 
u>(k) must be a 27r-periodic function of k (from now on 
the lattice constant is set 6x = 1). This can be incor- 
porated into the Lagrangian in Eq. ( |12| ) by the standard 
replacement d x (f> — > [<p(n + 5x) — 4>(n)[/8x, leading to the 
substitution k 2 — > 2[1 — cos(fc)] in the dispersion 



w(A)= 



(47) 



v/m 2 + 2a + 6b - (2a + 8b) cos(fc) + 26cos(2fc). 



The Green's function now reads 

dk 1 



„ ATTu(k) 



(48) 



The condition that the expression under the square root 
in u(k) is a complete square is again D(9) = with D 
defined in Eq. (|22| ) . When this is satisfied the dispersion 
simplifies to 



^vbs(fc) = myb B + 



^vbs 
TOvbs 



[1 — cos(fc)], 



(49) 



and l/w v bs(fc) has poles instead of branch cuts. [One 
pole within the Brillouin zone —n < Re(fc) < n with 
< Im(fc).] Now the second and fourth derivatives at 
k = are [cf. Eq. @] 



dk 2 {k 



m ok* 



^) = --,--■ (50) 



At the disorder point we find 
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=«)=-— 

dfc 4 m v bs 



(51) 



which is nonzero. The correlation length £ = l/Im(/c) at 
the disordered point is determined by the position of the 
pole, i.e., by the solution of the transcendental equation 



Ref. |T?j], Sf acting on the VBS ground state produces 
the linear combination of two states, one of which con- 
taining a boson at site i — 1/2, the other a boson at site 
i + 1/2. Hence the one-boson term of (SiSi+ n ) is in fact 

(Sffi) = 4(-)"[G (n - 1) + 2G^n) + G (n + 1)]. 



1 



Ovbs 



[1 - cos(fc)] = 0. 



(52) 



'vbs 



Using Eq. 



this leads directly to Eq. 



(55) 
including 



Working the other way around, knowing that at the VBS the (5-function piece, if 



(C-IC) point £ vbs = 1/ In 3, Eq. 
and thus the dispersion 



52) gives a vhs = 3m^ bs /2 



#0 = Wlvbs 



Wvbs = W v bs 



(53) 



With this expression the lattice Green's function defined 
in Eq. (ph reads 



G^n) 



1 



-\n\ In 3 



4m 



vbs 



(54) 



which should be conferred to Eq. (0). The functional 
form of the dispersion relation in Eq. ( p3| ) is exactly the 
same as the one appearing i.n the single mode approxima- 
tion of the VBS model.E3c3 There, one derives an upper 
bound for the gap A v b s = m-vbs = v40/9, whereas here 
we should use the phenomenological (numerical) value 
A vbs = 0.664314 in Eq. @. This, together with Eq. @ 
yields the value d 4 co/dk = — 3m v b s /2 ~ —1.0, which is 
rather close to the numerical estimate in Eq. (|4q). 

The split of the double root in the vicinity of the dis- 
order point can be analyzed similarly to the continuum 
theory. We do not go into details here, but emphasize 
that the critical exponents characterizing the behavior of 
and q(9) at # v bs, and the type of the singularity of 
the gap at #disp remain the same. Similarly, we find that 
the ground state energy is analytic everywhere. 

One can wonder about the possible consequences of 
keeping higher order terms in the continuum Lagrangian 
Eq. (fL2|) or in the improved lattice version. If the the- 
ory remains free the only effect is to bring about addi- 
tional branch cuts or poles in l/u){k). If the higher or- 
der terms are small the additional branching points are 
far in k space, and the C-IC transition remains intact. 
The long distance asymptotics of the correlation func- 
tions do not change. Since at the VBS point the exact 
correlation function in Eq. (0) only contains a single ex- 
ponential term, the free boson approximation does not 
allow any higher order spatial derivatives in the effective 
Lagrangian there. 

Although the dispersion w v b s in Eq. (|5^) gives rise to 
the exponential term in Eq. (Q), it misses the S n ^o con- 
tribution. This is, however, another artifact of the con- 
tinuum approach, which necessarily neglects some im- 
portant short distance details. In fact, we should recall 
that in a quantum spin liquid with short-range valence 
bond ground state the elementary excitations (bosons) 
are physically triplet bonds living between lattice sites, 
rather than on the sites themselves. As was argued in 



IV. SUMMARY AND DISCUSSION 

In summary, we proposed a simple effective field 
theory to describe the commensurate-incommensurate 
transition in the Haldane phase of the spin-1 bilinear- 
biquadratic chain. The theory is capable of reproducing 
many features of this transition previously seen in the 
numerical studies. Moreover it also has some new predic- 
tions. The effective theory predicts that the C-IC transi- 
tion at 9 v bs is not a phase transition in the conventional 
sense, since the ground state energy remains an analytic 
function of the control parameter 6. The only singularity 
occurring is in the correlation length. We should empha- 
size, however, that the correlation function itself remains 
analytic as a function of 9 for any fixed distance, unlike 
in conventional phase transitions. 

There is another point $disp close to the disorder point 
where another quantity becomes singular. This is the 
energy gap (Haldane gap) whose second derivative pro- 
duces a jump. This singularity in the singlet-triplet gap 
becomes important when a high enough magnetic field is 
applied, producing a crossing between these levels, and 
thus leading to the collapse of the gap. At the critical 
field, as 9 is varied a real phase transition takes place at 
#disp 7 thus this point is the endpoint of a phase transi- 
tion line on the magnetic field vs 9 plain separating two 
Luttinger liquid type phases. 

In a technical sense the C-IC transition is a conse- 
quence of an accidental degeneracy of roots of the disper- 
sion relation. We have shown in the free boson approach 
that this degeneracy causes the " dimensional reduction" 
of the correlation function, and makes it to be a pure 
exponential at the VBS point. We also found that the 
two-magnon contributions, and presumably any higher 
order, multi-magnon contributions too, vanish exactly at 
this point. We derived a formula valid in the vicinity of 
the VBS point showing how the pure exponential decay 
emerges from the standard form with algebraic prefac- 
tors. In particular we found that there is a crossover 
between a pure exponential decay and a decay contain- 
ing algebraic prefactors. The characteristic distance of 
this crossover tends to infinity as the VBS point is ap- 
proached. 

The spin-1 bilinear-biquadratic model studied in this 
paper is not the only model which produces a C-IC tran- 
sition. Another example is the spin-1/2 chain with next- 
nearest-neighbor interaction (this can also be visualized 
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as a two-leg zig-zag ladder) B By now it is well estab- 
lished that the Majumdar- Ghosh point of this model 
is a disorder point where a C-IC transition of the first 
kind occurs. On the same footing as described here, it 
seems possible to develop an effective theory which sup- 
poses that the elementary excitations (spin- 1/2 solitons 
in that case) are essentially free particles with a non- 
relativistic dispersion. Care should, however, be taken 
on the facts that solitons are always created in pairs and 
that they are spin-1/2 particles. Beside the S = 1/2 case, 
the appearence of disorder points has been demonstrated 
in other 5* > 1 frustrated Heisenberg chains tooEa An- 
other interesting quasi-one-dimcnsional system where a 
commensurate-incommensurate transition have been re- 
ported in a numerical investigation is-the SU(2)xSU(2) 
symmetric coupled spin-orbit model.E3 The elaboration 
and testing of effective theories, similar to the one de- 
scribed in this paper, for these models could be a possible 
direction of future research. 
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